function [ val ] = feature_dpratio(data, header, range, input_channels)

fftSize = 1024;
DC_margin = 10;

val = DP_ratio_measure_fn(data(range(1):range(2), input_channels), fftSize, DC_margin);

end

function [DP_ratio] = DP_ratio_measure_fn(data, fftsize, DC_margin)
% ratio of DC versus sidelobe location

DATA = fft(data, fftsize);
ML_peak = max([max(abs(DATA(1:DC_margin, :))) ; max(abs(DATA(end-DC_margin+1:end, :)))]);
SL_peak = max(abs(DATA(DC_margin:end-DC_margin+1, :)));
DP_ratio = ML_peak ./ SL_peak;
 
end